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ABSTRACT 

We analyze the stellar mass-to-light ratio (M/L) gradients in a large sample of lo- 
cal galaxies taken from the Sloan Digital Sky Survey, spanning a wide range of stellar 
masses and morphological types. As suggested by the well known relationship between 
mass-to-light {M/L) ratios and colors, we show that M/L gradients are strongly corre- 
lated with colour gradients, which we trace to the effects of age variations. Stellar M/L 
gradients generally follow patterns of variation with stellar mass and galaxy type that 
were previous found for colour and metallicty gradients. In late-type galaxies M/L 
gradients are negative, steepening with increasing mass. In early-type galaxies M/L 
gradients are shallower while presenting a two-fold trend: they decrease with mass up 
to a characteristic mass of M* ~ IQ^^-^Mq and increase at larger masses. We compare 
our findings with other analyses and discuss some implications for galaxy formation 
and for dark matter estimates. 
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1 INTRODUCTION 

Colour and absorption line gradients are markers of stel- 
lar population variations within galaxi es, providing impor- 
tant clues to galaxy evo l ution (e.g. ISpola.or et al.l 2009: 



201 ll: iRawle et al.ll20ld : 



2010al . herafter T-l-10, 



Spolaor et al 



Tortora et al 



Gonzalez-Per ez et al. 
2010; Tortor a et ajf 
2OIII ). Although spectral features are less ambiguous 
probe s of stellar populat i on pr o perties than colours are 
(e.g. ISpolaor et all |2009| . I2OI0I : iKuntschner et al.l I2OI0I : 
iRawle et al.l I2OI0I ) , there are only small samples available 
outside the nuclear regions of galaxies. Fortunately, studies 
based on colours turn out to provide, on average, similar age 
and metallicity estimates to spectral studies (e.g. T-l-10), al- 
lowing for much larger surveys of radially-extended stellar 
populations. 

In T-l-10 we analyzed colour gradients within an effec- 
tive radius (-Reff) for an exte nsive sample of low -redshift 
galaxies based on SDSS data (jBlanton et al.ll2005l ). includ- 
ing both early-type (ETGs) and late-type (LTGs) galajdes. 
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We carried out simple stellar populations analyses, fitting 
synthetic models to the observed colours to derive estimated 
age and metallicity gradients, and stellar masses (see T-f 10 
and later in the paper for further details). The main results 
found can be summarized as follows: 



(i) Colour gradients are, on average, negative (i.e. galax- 
ies are redder in the centre), except for galaxies with lower 
stellar masses (M* < 10^'^ Mq) which tend to have null or 
positive gradients. At fixed mass, LTGs have steeper gradi- 
ents than ETGs, whose colour profiles in general show milder 
variations . 

(ii) The colour gradients of LTGs decrease systematically 
with mass while the trend for ETGs changes near a mass of 
M« ~ 10^"'^ Mq: systematically decreasing and increasing 
at lower and higher masses, respectively. 

(iii) The primary driver of the colour gradients, mainly 
for galaxies with older centers, is metallicity variations, 
with age gradients playing a secondary role. For LTGs, the 
metallicity gradients tend to be steeper for higher masses. 
The same is true for ETGs except that for masses above 
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Figure 1. Diagram of colour versus stellar mass for the local 
SDSS galaxy sample analyzed in T+10. The grey-scale contours 
show the density of data points, with brighter regions marking 
higher densities. The thin white line separates the red sequence 
and blue cloud. The arrows and the curves in the inset panel give 
information about the efficiency of the phenomena which drive 
the two fold trend we discuss in the paper (see T+IO for further 
details). AGN feedback and merging are important at high mass 
with an efficiency that increases with mass, while SN feedback 
and gas inflow drive the galaxy evolution at lower masses, with 
an efficiency that is larger at the lowest masses. 



M* r^ 10^"'' Mq, the (negative) gradients are fairly invari- 
ant or slightly increasing with mass. 

(iv) At fixed mass, galaxies with older centres have the 
shallowest metallicity and age gradients. 



We found that these results we re generally consis- 
tent with others i n the literature (e.g. ISpolaor et al.l l2009l : 
iRawle et al.l I2OI0I ) once different galaxy selection effects 
were taken into account. 

The overall observational picture can be framed within 
the predictions of hydrodynamical and chemodynamical 
simulations o f galaxy formation IIBekki fc Shioval 199£ ; 
Kawatal I 2OOI": Kawata & Gibson 2003; Kobayashi 200j ; 
Hopkins et al. 2009 : Pipino et al. 2010 : Tortora et al. 2011 ) 



which include various physical processes that can influence 
the gradients: supernovae and AGN feedback, galaxy merg- 
ers, interactions with environment, etc. A fresh and quali- 
tative overview of this physical scenario is given by Fig. [T] 
where we show the isodensity contours for the colour-mass 
diagram (which we have already discussed in T-f-lO) and 
superimpose some lines (in the inset) and arrows which in- 
dicate the efficiency of the highlighted physical processes. 
In particular, the observations for high-mass ETGs seem to 
generally suggest t he influence of merge r s and AGN feedbac k 
l|Kobavashil I2OO4I : iHopkins et all [2OO9I : ISiiacki et al1l200'i1 ). 
with an efficiency possibly increasing with mass (see red 
arrows and lines in the inset panel). On the other hand, 
the population gradients of LTGs and less massive ETGs 
may be mainly driven by infall a nd supernov ae feedback 



llLarsonlll974 (l975l: lKawatall200ll: iKawata fc G ibson 200p; 
iKobavashil |2004| : iPipino et all I2OI0I : iTortora et all mm . 



which are more efficient at very low masses (see blue lines 
and arrows). 

In this paper we consider the implications of the ob- 
served colour, metallicity and age gradients from T+10 for 
variations of stellar mass-to-Ught ratio M/L with radius. 
This issue is relevant for understanding the detailed dis- 
tribution of mass within galaxies: both the stellar mass 
and the dark matter, whose distribution is inferred by sub- 
stracting the stellar mass from the total dynamical mass. 
There has been scarcely any literature to date that fo- 
cu ssed on this topic, e xcept for a study in spiral galaxies 
bv lPortinari fc Saluccil ([2OIO). 

As a first approximation, one could u se single-colour 
semi-e mpiri cal model s such as those of iBell fc de Jona 
(|200ll) and iBeU et al.l l|2003l . BOS hereafter) to go directly 
from colour gradients to M/L gradients. However, we wish 
to make use of multiple colours, to allow for departures of 
individual galaxies from mean trends, to explore the sys- 
tematics thoroughly, and to trace the dependence on stellar 
parameters, like age and metallicity. 

This paper is organized as follows. The sample and the 
spectral analysis are discussed in ij2] in where we also show 
the ability of our fitting procedure to recover the correct 
M/L gradients, and we examine possible systematics. The 
results are shown in ^ while 31]is devoted to a discussion 
of the physical implications and conclusions. 



2 DATA AND SPECTRAL ANALYSIS 

Our database consists of 50 000 low redshift (0.0033 ^ z ^ 
0.05) galaxies in the NYU Value-Added Galax y Catalog 
extracted from SDSS DR4 (|Blanton et al.l |2005| , hereafter 
B05i3, recently analyzed in T-|-lcB Following T-l-10, we 
have identified ETGs as those systems with a r-band Sersic 
index satisfying the condition 2.5 ^ n ^ 5.5 and with a con- 
centration index C > 2.6. The final ETG sample consists 
of 10 508 galaxies. Of the remaining entries of our database, 
- 27 800 are LTGs, defined as objects with C sj 2.6, n sj 2.5, 
and (To ^ 150 km/s. 

We have used the multi-band structural parameters 
given by B05 to derive the colour profile {X — Y){R) of each 
galaxy as the differences between the (logarithmic) surface 
brightness measurements in the two bands, X and Y . We 
have defined the colour gradient as the angular coefficient of 

c/ Y" V\ 

the relation X - Y vs \ogR/RoH, Vx-y = t^ — 77777^ — T' 

S\og{R/Rcitj 

measured in mag dex~^ (omitted in the following), where 
Raff is the r-band effective radius. By definition, a positive 
gradient, Vx-y > 0, means that a galaxy is redder as R 
increases, while it is bluer outward for a negative gradient. 
The fit of synthetic colours is performed on the colours at 
-Ri = -Rcft/10 and R2 = -Reft and on the total integrated 
colours. 

The stellar population properties (age, metallicit y and 
M/L) are derived by fitting lBruzual fc Chariot I l|2003l , here- 
after BC03) "single burst" synthetic stellar models to the 



^ The low redshift NYU-VAG catalogue is available at: 
http : //sdss . physics . nyu . edu/vagc/lowz . html 
'' See T-(-10 for further details about sample selection, incom- 
pleteness and biases. 
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Table 1. Recovery of stellar population parameters from Monte Carlo simulations. A purely random sample of galaxies is generated, 
leaving the ages and metallicities free to vary across the whole range. We list the median error and the 25-75 percent scatter in the 
gradient A(X) = Xfjt — X^^, where X^^ and Xflt are the input and fitted gradients with X = Vt», Vagc and Vz. The data are divided 
in three subsamples, on the basis of the input M/ Ls at i?2- Different initial perturbations 5 and the results (1) with optical only (ugriz) 
(2) optical and near-IR (ugrizJHKs), and (3) optical and UV (ugrizFUV NUV) are shown. 
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Table 2. Median errors as in Table [T] but for 5 ■■ 
values. 
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observ ed colours . A^e and metallicity are free to vary, 
and a IChabrier I (|200in initial mass function (IMF) is as- 
sumed. Stellar parameter gradients are defined as Vvi^' = 



where W 



are the estimated 



6\o^{R/R,sy 

age, metallicity and M/L. Because of the definitions adopted 
for Ri and R2, an easier and equivalent definition can 
be used, in fact, we define Vagc = log[age2/agei], Vz = 
Iog[Z2/.Zi] and Vx. = log[T«2/T*i] where (age,, Z,, T.^), 
with i = 1,2, are the estimated parameters at R\ and 
7?2, respectiveljo. In the present paper we will concentrate 
on the analysis of the M/L gradients, Vt, , in terms of 
stellar mass and velocity dispersion. The dataset, includ- 
ing main stellar parameters and gradients, is available at: 
[http : //www ■ itp . uzh . ch/~ctortora/gradient_data . html 
Before proceeding, we check for systematic effects on 
our Vt, estimates from the stellar populati on fits. Because 
of the well known age-metallicity degeneracy |Worthevl 19941 . 
iBruzualfc Chariot! 12001 iGallazzi et all l2005t ). the stel- 
lar population parameters and the gradients might be bi- 



We use B-band M/L, but for simplicity we omit the subscript 



ased when using the optical colours only, as we do now. 
Widening the wavelength baseline to include near-infrared 
(NIR) colours should ameliorate the degeneracy. In T-l-10 we 
checked age and metallicity inferences using optical versus 
optical-|-NIR constraints, and found little systematic differ- 
ence. Here we will carry out a similar analysis of the M/L 
gradients, using a suite of Monte Carlo simulations. 

We extract some sets of simulated galaxy spectra from 
our BC03 spectral energy distribution libraries with random, 
uncorrelated stellar parameters and gradients. We start from 
generating data with fully random ages and metallicities, 
and then we also analyze the cases when the stellar pa- 
rameters are constrained within fixed intervals, generating 
not-zero sample average gradients. We extract the input 
galaxy colours and add simulated measurement errors, as 
randomly extracted steps from the interval (— 5,-|-(5), with 
5 — 0.01, 0.03, 0.05. We apply our fitting procedure, search- 
ing for the best model in our synthetic library, which re- 
produces the colours of each of the simulated galaxies and 
then compare the output parameter estimates to the input 
model values. We perform the fit (1) using only the optical 
SDSS b ands ugriz and (2) adding NIR photometry (J, H 
and Jf5: |jarrett et al.ll2003l ) and (3) ultraviolet photometry 
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(NUV and FUV: iMartin et al.ll2005l ). We define the median 
error in the gradients as A(X) = Xfit — ^in, where Xin and 
Xflt are the input and fitted gradients with X = Vt» , Vagc 
and Vz, which we show in Tables [T] and [2] 

Starting fi:om the purely random case shown in Table [T] 
we have found that, adding the NIR or UV data, both the 
M/ L values and their gradients are perfectly recovered with 
a very little scatter. With optical data only, the uncertain- 
ties are doubled (around 0.1 in A(Vt.)) and there is a small 
systematic shift in the gradients (of log Vt. ^ 0.01 — 0.02 for 
the worst case analyzed), but the results are fairly well recov- 
ered with no spurious trends expected. We notice that the 
median results are quite independent of the range adopted 
for the input T,2, and only slight differences in the scat- 
ter are found. Although we will be mainly interested in the 
trends of Vt, with mass and velocity dispersion, we have 
also checked the systematics in age and metallicity gradients. 
We find similar results to the case of Vx, discussed above, 
but slightly larger shifts in the median and larger scatters 
(~ 0.2 in the worst cases) when only the optical is used. In 
Table [2l calculating the sample averages, we start with the 
input values Vago 7^ and Vz 7^ 0. Here the median er- 
rors and the scatters increase, but the conclusions discussed 
above are qualitatively confirmed. When optical data are 
used, for high T*2 we have the worst discrepancies in the 
estimated M/L gradients amounting to ~ —0.05, —0.07. For 
age and metallicity gradients, we could underestimate (over- 
estimate) the age (metallicity) gradients of ~ 0.1, which is 
reduced to ~ 0.02 when near-IR or UV are included. 

Thus, in many cases we have found that the median 
errors are acceptable and within the typical sample scatter 
and the uncertainties due to the adopted stellar population 
prescription. We will consider these results as qualitative 
indications of the effects on estimated gradients caused by 
systematics. 



3 RESULTS 

In this section we discuss first the M/L gradients in terms of 
colour, age and metallicity gradients and then as a function 
of stellar mass and velocity dispersion, paying attention to 
the role of central age in the observed correlations. 



3.1 Dependence on colour, age and metallicity 
gradients 

We start by plotting in Fig. [2] the M/L gradients as a 
function oi g ~ r colour gradient (Vg_r) for ETGs and 
LTGs. The two quantities are strongly and positively cor- 
related and no statistically significant difference between 
the morphological types is found. While such a correlation 
is not surprising owing t o the well known-cor relations be- 
tween colours and M/L (|Bell fc de Jong||200ll : B03), which 
should apply to all radii thus producing the relationship be- 
tween Vt. and Vg-r , there are important details to con- 
sider. iBell fc de Jong (l200l|) have adopted different evo- 
lution models for LTGs to establish some linear correla- 
tions between integrated total colours and the corresponding 
M/L determining some best fit relations between these two 
quantities, which have been proven to work for ETGs as 



well IIBell et al.l 20031) . Adop ting the solar metallicity mod- 
els in iBell fc de Jong 112001^ w e have fou nd a good agree- 
ment with iBruzual fc Charloti l|2003l ') or iMarastonI (|2005| ) 
synthetic prescriptions in ETG s with high M/L, while at 
low M/Z/ iBell fc de Jond l|200ll ') overestimates the M/L (see 
iTortora et al.l ([20091) for further details). In the left panel 
of Fig. [2] we compare our VT,-Vg_r trends with the ones 
(solid and dashed lines) derived from the BOS best fitted re- 
lations, finding a good agreement. The black and green lines 
plotted in Fig. [2] are derived by the relations between the 
g-band M/L and the colours (u — g, g — r, g — i and g — z) 
from BOS. We insert in these relations the measured colour 
for each galaxy in our sample, deriving the corresponding 
M/L. This operation is made for both the radii _Ri and R2, 
thus we can calculate the M/L gradients. Finally, the plot- 
ted lines are obtained, first calculating medians of Vx, in 
bins of Vg_r and then fitting a straight line. We see that the 
best fitted relations depend on the colours adopted in the 
colour-M/L relations from B-I-OS to derive the M/L gradi- 
ents. For instance, we observe a fiattening of the Vx,-Vg_r 
relation when using colours probing redder wavelengths, i.e. 
g — i and g — z, but a shallower trend is obtained when u — g 
is adopted. The agreement is quite good for moderate colour 
gradients, while for very steep positive or negative gradients 
some departures are observed. 

Despi te the good ag r eemen t, the very simplified ap- 
proach in iBell fc de Jong (120011 ) and BOS has some strong 
limitations. In fact, a) the BOS relations hold for integrated 
quantities, thus the radial variation of the M/L can be de- 
rived if one assumes that the same relations apply at differ- 
ent apertures and in general, the zero-point and the slope 
of these fitted linear relations would depend on the radius 
adopted and b) these relations do not take into account 
the stellar parameters (e.g. age and metallicity at each ra- 
dius) thus they do not allow to catch the full physics behind 
these correlations. These shortcomings are overcome by our 
multi-colour direct stellar population analysis which con- 
nects the M/L gradients self-consistently to the underlying 
stellar population properties (such as age and metallicity). 
This also allows us to track the physical reasons for the M/L 
variations, which in turn should be reproduced by galaxy 
formation models. 

Differently from what happens for colour gradients 
which are mostly driven from metallicity gradients, from the 
right column of Fig.[2]it is evident that the Vx.-Vg-r rela- 
tion is driven mainly by age, in both LTGs and ETGs. 



3.2 Trends with stellar mass and velocity 
dispersion 

The trends of Vx, with stellar mass and velocity dispersion 
are shown in the top panels of Fig. |3] As with the trend of 
colour and metallicity gradients with mass (see comments 
in !JT] also T-j-lO), we find that LTGs have a decreasing 
trend with mass and show Vx, > for log M^/Mq < 9 and 
Vt,< at larger masses. In contrast, ETGs have positive 
Vx, ~ 0.25 at low masses, which decreases down to a min- 
imum of Vx, ~ —0.1 at logM,/MQ ~ 10. S and then invert 
the trend, with the most massive ETGs having gradients 
around zero or slight positive. When plotted as a function of 
velocity dispersion, the Vx, run with ctq is generally fiatter, 
with LTGs having no trend and ETGs showing a two-fold 



M/L gradients 5 



0.5 



-0.5 



i i 



-0.6 -0.4 -0.2 

Vg-r 



0.2 



0.4 



0.6 




Figure 2. Left. M/L gradient as a function oi g — r gradient for ETGs (red squares) and LTGs (blue triangles). Medians of the gradients 
(in bins of colour gradient) are plotted together with the 25-75 per cent quantiles, shown as error bars. The lines are for the gradients 
obtained using the best fitted M/L-colour relations in BOS. The relationships between the g-band M/L (g-band is approximatively the 
same as the B-band) and g — r (solid black line), u — g (long dashed green line), g — i (short dashed green line) and g — z (green line) are 
adopted and the recovered Vt, are plotted as a function of Vg_r. Right. Age (top) and metallicity (bottom) gradients as a function of 
g — r colour gradient. Solid and dashed lines show systems whose centres are older and younger than 6 Gyr, respectively. M/L and colour 
gradients are tightly correlated for both ETGs and LTGs, and, while Vt» 
Vz. 



trend with a minimum around logao = 2.1 km/s and a mild 
increase toward the low and high cro sides, always staying 
consistent with Vt, < 0. 

Since galaxies are morphologically selected by using cri- 
teria on concentration and Sersic index, the steeper gradi- 
ents found for LTGs when compared with ETGs indicate a 
sequence in terms of such structural parameters. In fact, at 
fixed mass, galaxies with lower C or n have steeper colour 
and M/L gradients. 

As observed for the stellar population gradients (see 
T-I-IQ), the central age is responsible for much of the scatter 
of the trends versus mass and velocity dispersion, which is 
illustrated by the median distribution for centrally old (> 
6 Gyr) and young (< 6 Gyr) galaxies in Fig. |3l In particular, 
young LTGs have almost null gradients, while the oldest ones 
have very steep values of Vt. ~ —0.3 to —0.4. ETGs show a 
similar but oflset effect. Older ETGs have slightly negative 
values of ~ —0.1, while the younger ones have steep positive 
Vt, ~ 0.5 at very low masses, decreasing to ~ 0.1 at larger 
masses. 

As a comparison, we also show in Fig. [3] (top panels) the 
results when galaxy colours are fitted to a synthetic spectral 
model which has no age gradient (fixed age = 10 Gyr). In 
this case, LTGs have Vt, which are null up to log Mt/Mg ~ 
9.5, and decrease at larger masses, while gradients in ETGs 
show a milder trend when compared with the reference fit. 
When plotted versus ctq both LTGs and ETGs have similar 
Vt, ~ —0.1 for all velocity dispersions. Note that these 
results are consistent with what was found for the older 
galaxies in the reference fit. We also compare our results 



with g-band M/Ls derived from BOS. As an example we 
plot the median lines derived using the AI/L — (g — r) and 
M/L - [u - g) relations for both ETGs and LTGs. The 
trends with M, and velocity dispersion are qualitatively un- 
changed, but some differences emerge since the M/Ls from 
each BOS relation are derived using different colours, which 
probe different wavelength spectral regions. The M/Ls ob- 
tained using the g — r colours seems to resemble our trends 
for LTGs, while the agreement is poorer when u — g is used. 
For ETGs our trend is between the two BOS results at low 
mass, while the M/L from u — g is a. little bit better at large 
mass. 

We have also checked the effect of the waveband 
adopted for the M/L ratios in the bottom panels of Fig. [S] 
Here, the B-band M/L gradients are compared to the ones 
in V- and i-band. Except for the most massive ETGs, in 
general the M/L gradients tend to become shallower (closer 
to zero) when using redder bands. Also, the trends with 
mass and velocity dispersion become weaker. Thus, redder 
bands, which are less affected by dust extinction and trace 
the older and metal-richer stellar populations in galaxies, 
tend to be less dependent on the existence of gradients in 
colour and stellar population parameters. Extending such 
analyzing into the NIR bands should find M/L gradients 
increasingly approaching zercQ. 

From the preceding results, we can see that although the 



We prefer to limit the analysis to the wavelength coverage 
available for our dataset, avoiding extrapolations to near-IR wave- 
bands. 
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Figure 3. M/L gradients as a function of stellar mass (left) and velocity dispersion (right). For ETGs a two fold trend in terms 
of both M* and ctq is found. LTGs present a steepening of the gradients with M«, while they are independent of ctq. Top: Default 
B-band gradients. The symbols are as in Fig. [2] with solid and dashed lines indicating galaxies with central age older and younger than 
6 Gyr, respectively. The light and dark grey lines are the results for ETGs and LTGs when no age gradients are assumed (adopting 
age = 10 Gyr). Middle: Default _B-band gradients, symbols are as panels above. Black and green lines are the medians obtained using the 
M/L — {g — r) and M/L — [u — g) relations from BOS, continue and dashed are for ETGs and LTGs. Bottom: Average gradients using 
different bandpasses: B (red and blue lines for ETGs and LTGs, respectively), V (orange and violet) and i (yellow and cyan). 



age gradient is generally the driving parameter behind the 
M/L gradient, this conclusion can vary depending on galaxy 
type, mass and central age. We summarize the gradients for 
different galaxy subsamples in Table [S] noting also which of 
the stellar parameters drives the M/L gradient. 



4 DISCUSSIONS AND CONCLUSIONS 

In this paper we have investigated the correlation between 
M/L gradients and stellar mass or velocity dispersion for a 
large sample of local galaxies from SDSS. The gradients Vx, 



have been derived by the fitting of synthetic spectral models 
to observed optical colours. We have used Monte Carlo sim- 
ulations to check that the optical-only wavelength baseline 
should not dramatically affect the results compared to in- 
cluding near-infrared and ultraviolet data. This conclusion 
is encouraging in view of future studies of higher redshift 
galaxies which are generally observed in rest-frame optical 
bands. Anyway, future analysis in local samples, relying on 
near-IR and/or UV data would be able to give further de- 
tails about the degeneracies and systematics affecting stellar 
fitting. 
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Table 3. M/L, age and metallicity gradients for subpopulations of galaxy type, stellar mass and central age. In (parentheses) we note 
whether the age or metallicity gradient drives the M/L gradient. 



LTGs 


logM./Mg <9 




logA/./Mo ^9 




Young 


Vt. > 0; Vagc > 0; 
(age) 


Vz <0 


Vt. < 0; Vagc > 0; 
(metallicity) 


Vz <0 


Old 


Vt. < 0; Vago < 0; 
(age) 


Vz<0 


Vt. < 0; Vagc < 0; 
(metallicity) 


Vz <0 


ETGs 


logM*/MQ < 10 




logA/./A/Q ;>10 




Young 


Vt. > 0; Vagc > 0; 
(age) 


Vz;>o 


Vt. > 0; Vagc > 0; 
(age) 


Vz <0 


Old 


Vt. <0; Vago < 0; 
(age) 


Vz^O 


Vt. < 0; Vagc ~ 0; 
(metallicity) 


Vz <0 



We have found that there exists a tight positive corre- 
lation between Vt. and colour gradients which is not sur- 
prisi ng because of the wel l known Af/L-colour correlations 
(e.g.. lBell fc de Jondliooll : BOS). We have found that M/L 
is a decreasing function of stellar mass in LTGs, while a 
two-fold trend is found for ETGs (for both stellar mass and 
velocity dispersion), in agreement with the trends of colour 
gradients discussed in T-l-10. Similar results are obtained 
using the M/L-colour relations from B03, which, although 
easy to use, do not give information about stellar popula- 
tion parameters. Thus, the novelty of our approach resides 
in our ability to check how Vt. gradients are driven by vari- 
ations in metallicity or age. In general, we have found that 
age is most important, but that metallicity becomes more 
important in many high-mass galaxies. 

This picture does not change dramatically if M/L gradi- 
ents are calculated in redder bands {V- and i-bands rather 
than B). In these cases Vt, are shallower and the trends 
with stellar mass and velocity dispersion become milder 
(Fig. O. These results suggest that M/L estimates from 
red or infrared bands would not depend on galactocentric 
radius, since the radiation emitted by stellar populations in 
these spectral regions is more homogenous. 

This is the first work to derive M/L gradients 1) from 
observations of a large galaxy sample and 2) using a di- 
rect stellar population analy sis approach. The only po int 
of comparison is the work bv lPortinari fc Saluccil ([2010) in 
a sample of spiral galaxies, where negative M/L gradients 
were found to be important in modelling rotation curves. 

We may explain some of our observed trends through 
physical mechanisms as follow. Starting from low-mass sys- 
tems, centrally-young LTGs have age increasing with radius 
while metallicity decreases, which can be due to the presence 
of a gas inflow which allows the formation of younger stars in 
the centre with a smaller M/L. Centrally-old LTGs have age 
decreasing with radius which implies recent outer star for- 
mation with a similar metallicity to the centre. Here the cen- 
tral star formation might have been quenched by SNe which 
have expelled metals to larger radii (where the younger stars 
drive a smaller M/L). Centrally- young ETGs have positive 
age, metallicit y and M/L grad ients, and may hosting ex- 
panding shells (|Mori et al.lll997l 'l. Centrally-old ETGs show 
mild negative age gradients and very shallow positive metal- 



licity gradients, similarly to the old LTGs thus again sug- 
gesting SNe having played a role. 

The more massive centrally-young LTGs have decreas- 
ing M/L with radius, but also strong negative Vz, thus 
the new stars formed inside the core (Vage> 0) have higher 
metallicity and M/L. This might be due to recycled gas 
from SNe or AGN which is falls back to the centre due to 
the deep potential wells and produces a rejuvenated stellar 
population. In centrally-young ETGs, on the other hand, 
the younger stars formed later in the centres (i.e. with 
positive age gradients) and with higher metallicities, pos- 
sibly becaus e of a central burst of star formation in "wet" 
mergers (e.g.lRothberg fc Josephll2004l : iHopkins et"alll2009l : 
iHopkins et al.ll201(]| : T-HOV Centrallv-old LTGs and ETGs 
show very similar behaviors, which seem compatible with 
merging events that largely wash out any age gradients (and 
the Vz in ETGs), and may produce coeval metallicities in 
the outer regions of the LTGs. In both cases the net effect 
is that metallicity is the main driver of the strong negative 
Vt, in LTGs and of the weak negative Vt. of the massive 
ETGs (in any). Thus, the physical mechanisms that can ac- 
count for the observations outlined in Tab.[3]seem consistent 
with the evolutionary scenario discussed in SJT] 

In particular, the two-fold trend of ETGs in Fig. 
[3l with a minimum around XogM^/ Mq = 10.3 is mir- 
rored by the colour and stellar population gradients as 
in Fig. [T] (and T-j-lO) and found a l so in other works 
(JRawle et al.ll2010l : [Spolaor et all [20091 . \20ldi ). The charac- 
teristic mass agrees well with the t ypical mass scale break fo r 
star forming and passive systems l|Kauffmann et aini2003l'). 
for "b r ight" and "ordinar y " gal a xies (ICapaccioli et al.l 
'1992a; Graham fc Guzman 2003; 'Grah am et all l2003l : 
Truiillo et al. 2004 : Kormendv et al. 2009). or for cold and 



hot flow (shock heating: iDekel fc Birnboimll2006l ). 

Following the main conclusions in T-l-10, quasi- 
monolithic collapse and supernovae feedback could be 
the main drivers of the de creasing Vt. - M. trends in 
LTGs and low m ass ETGs (JLarsonl [l97i : iKawatal 1200 ll . 
iKawata fc GibsonI [2003), while the increasing trends in 
massive ETGs could be explained by th e heightened im- 
portance of mergers and AGN feedba ck (|Kobavashil 12004 : 
ISiiacki et a"i]|2007l : [Hopkins et aLlbOOgD . 

M/L gradients are expected from simulations of galaxy 
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formation and evolution. E.g. iBoissier fc Prantzos I (| 19991 ) 
found in chemo-dynamical simulations of spiral galaxies that 
the central regions (< 12 kpc) have steep negative M/L gra- 
dients, becoming flatter at large r radii. High-redshif t {z ~ 2) 
merger remnants simulated by IWuvts et all (|2O10l ) showed 
time-dependent negative M/L as a consequence of redder 
cores in redder galaxies (similarly to observations of local 
galax;ies in T-(-10). M/L gradients were also predicted for 
merger remnants bv lHopkins et al.l l|20ld ). who stressed that 
the effective radius derived from the _B-band light would in- 
crease with time (by ^ 10%) simply because of stellar pop- 
ulation gradients. This effect is important for conclusions 
about stellar densities and related dynamical implications 

i e.g. for the fundamental plane). Similarly Ivan DokkumI 
200a) argued that the existence of such gradients can imply 
that the observed quiescent galaxies at large redshift may be 
of an order or magnitude more dense than previously found 
from observations. 

The M/L gradients present a general complication 
for estimates of th e dark matter distribution in galaxies. 
IPortinari fc Saluccil (|201G| ) found that for spiral galaxies, in- 
cluding M/L gradients would change the fits to circular ve- 
locity curves but not produce strong effects on the inferred 
dark matter density parameters. To our knowledge, this is- 
sue has not been addressed in the literature for early-type 
galaxies, and it is beyond the scope of the present paper to 
do so in any detail, fn particular, we have studied the stellar 
population gradients only inside 1 -Rcff, while the gradients 
out to ~ 5 -Reff could b e importan t for d ark matter analyses 
fe.g. iNapolitano et al.i r2005. 200d l201ll 'l. Assuming that the 
gradients do not change radically outside ~ 1 RcS, for the 
massive ETGs generally studied with lensing or dynamics 
((To f> 150 km s~^), we would expect the stellar M/L gra- 
dients to be mild (see Fig.[3|, and to not strongly affect the 
dark matter estimates. 

A final issue is that until now, we have assumed that 
the IMF is invariant. If there are galaxy-to-galaxy var i ations 
in the IMF (see e.g. Ivaii Dokkum fc Conrovl I2OI0I . I2OIII : 



iGun awardha na et al.ll2011l ). then the stellar M/L values are 
changed only by a constant multiplicative value, and the 
implied gradients remain unchanged. However, IMF varia- 
tions if real would then likely exist within galaxies owing 
to variations in stellar populations. Indeed, there have been 
many suggestions that stjirs forming at early times had a 
non-standard IMF (e.g. jLarson 20051: iKlessen et all l2007l : 
Ivan Dokkumll2008l : [Pat^boOSl : iHolden et al.ll2010l ) . 

One study of nearby ETGs suggested that systems with 
older stars within 1 R^s have "lighter" (closer to Ghabrier 
than Salpeter) IMFs C Napolitano. Romanowskv fc Tortora I 
I2OIO1 V This would imply that positive and negative age gra- 
dients would decrease and increase the stellar M/L gradi- 
ents, respectively. Revisiting Fig. O the gradients for the 
centrally-young objects would decrease and the centrally old 
objects would increase, possibly reducing and even eliminat- 
ing the age-related scatter in the gradient trends. 

Another scenario would be if the central regions of 
massive ETGs have "higher -mass" IMFs than standard 
l|van Dokkum fc Gonrovll2011l ). while the outer regi ons have 
accre te d galaxies with more "normal" IMFs (e.g. iKroupal 
I2OOII or IChabrieiil2003l ) . This would decrease these galaxies' 
stellar M/L gradients, although it is not clear how this effect 
would relate to age and metallicity trends. 



In future analyses, we plan to analyze stellar M/L pro- 
files out to a few i?cff, and analyze the impact of M/L gra- 
dients on dark matter p rofile inferences, wh ich may be com- 
pared to AGDM (e.g. iNavarro. Frenk fc White 1996). We 
will also consider the (probably small ) adjustments implie d 
for central dark matter inferences (e.g. lTortora et al.ll2010a ). 
Dust effects will further be examined: dust is not important 
for metallicity gradients (see T+10) but could affect age and 
thereby M/L gradients, particularly for intermediate-mass 
ETGs. 
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